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Abstract 

We show that the quark fragmentation function D(z) and the quark distribution 
function q(x) are connected in the z — > 1 limit by the approximate relation D(z)/z ~ 
q(2 — 1/z), where both quantities are in their physical regions. Predictions for proton 
production in inelastic e + e~ annihilation, based on the new relation and standard 
parametrizations of quark distribution functions, are found to be compatible with 
the data. 
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Inclusive deep inelastic scattering (DIS) and hadron production in e + e~ inelastic 
annihilation (IA) are important sources of information on the structure of hadrons. 
By simple crossing it is possible to relate the structure functions of these two processes 
[ED]. However, the relation so obtained (known as the Drell-Levy-Yan relation) is of 
little use, as it connects the structure functions in the physical region of one process 
to the structure functions in the unphysical region of the other process. 

It would obviously be important to have a relation between the DIS and IA struc- 
ture functions, taken both in their physical regions. We could then exploit the accu- 
rate information we already possess on the quark distribution functions of the nucleon 
to predict the quark fragmentation functions, which are still poorly known, or vice 
versa, by measuring the fragmentation functions of hadrons which cannot be used as 
DIS targets we could predict the quark densities inside those hadrons. An example 
which immediately comes to mind is the A: the quark dynamics inside this hyperon 
is highly relevant for our understanding of the spin and flavor structure of hadrons 

a as- 

As a matter of fact, a relation connecting the structure functions of DIS and I A in 
their physical regions does exist in the literature. It is the so-called Gribov-Lipatov 
"reciprocity" relation. As we shall see, this relation, in its commonly used form, has 
no real justification. Moreover it is not supported by phenomenological evidence. The 
purpose of this paper is to derive another, well founded, relation connecting DIS and 
IA, and to show that this relation, within its range of validity, is in good agreement 
with the existing data. 

The DIS cross section is written in terms of two structure functions F\(x, Q 2 ) and 
F 2 (x, Q 2 ), where x = Q 2 /2p-q is the Bjorken variable and Q 2 = —q 2 is the momentum 
transfer squared. Two analogous quantities appear in the IA cross section: they are 
denoted by Fi(z, Q 2 ) and F 2 (z,Q 2 ), where now z = 2p ■ q/Q 2 and Q 2 = q 2 is the 
center-of-mass energy squared. In leading order QCD the Callan-Gross relations 
connect F\ (Fi) to F 2 (F 2 ) as 

F 2 (x,Q 2 ) = 2xF 1 (x,Q 2 ) } (1) 
2zF 1 (z,Q 2 ) = -z 2 F 2 (z,Q 2 ), (2) 

and, at the same order, the structure functions can be expressed in terms of the 
distribution functions q a ,a and the fragmentation functions D a ^ as 

2F 1 (x, Q 2 ) = -F 2 (x, Q 2 ) =J2el Q 2 ) + Qa(x, Q 2 )} , (3) 
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2zF 1 {z, Q 2 ) = -z 2 F 2 (z, Q 2 ) = 3 ]T e 2 a [D a (z, Q 2 ) + D^(z, Q 2 )] , (4) 

a 

where the sums run over all flavors (the factor 3 comes from a sum over colors). 
The traditional form of the Gribov-Lipatov relation reads || 

zF l {z) = F l {z), 
z 3 F 2 (z) = -F 2 (z), 

where F\ )2 (z) means that the DIS structure functions are evaluated at x — z. Phe- 
nomenological tests of this relation have been carried out || [7|, || and it turns out 
that the IA structure functions predicted by Eq. (|5]) undershoot the data. But the 
main shortcoming of (^) is its uncertain theoretical status. In fact, what Gribov and 
Lipatov proved in their classical papers |9| is that the non-singlet splitting functions 
for DIS and IA are equal at leading order (for a detailed and clear discussion see [ |i~0| ) 

T l £(z) = P£\z). (6) 

Thus Eq. (|5|) is true at leading order if one assumes that at a nonperturbative scale 
/i 2 the input distribution and fragmentation functions are related to each other by 

3D(z) = q(z) . (7) 

This relation is unjustified. We will now show that the true nonperturbative relation 
existing between q(x) and D(z) in the large- z limit is the approximate relation 

-D(z) ~ q(2-l/z) . (8) 

z 

Since 2 — 1/z ~ i-(i-\/ z ) = z, as z —* 1, Eq. (g) can be approximated further by 



-D(z) ~ q(z) . (9) 

z 



This relation was used as a phenomenological Ansatz in M. 1 ] 



Let us start from the general definition of the quark distribution function [12 



q(x) = J ^e^-(h(p)\mi + ^(r)\h( P )), (10) 

for a hadron h with mass M and momentum p = (p°,p). The light-cone compo- 
nents are defined as p^ = \{p° ±p 3 ). The normalization of the states is (p\p') = 
(2ir) 3 2p + 5(p + — p' + ) S 2 (p ± — p'jj. By inserting a complete set of intermediate states 
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|n) with momentum p n = (p^,p n ) and mass M n and making use of the translational 
invariance, one obtains 

^) = ^E / ^5(pl-M*)5(p+-xp+- P mn(Pn)\Mm(p))\ 2 , (11) 

where ip + = §7~7 +, is the 'good' component the quark field operator. 

Similarly, the fragmentation function of a quark into an unpolarized hadron h is 
defined as |12j (a sum over the final spin of the hadron is performed) 

X Tr{7+(0|^(0)|Mp),n(p„))(Mp),«(p„)W'(r)|0}, (12) 

where the initial quark carries a light-cone momentum k + = p + / z. The normalization 
of D(z) is such that J2h I dzzD(z) = 1. Using translational invariance one obtains 


-D(z) = -1=5: / ^^Pl-M 2 J6( P yz- P + -p^\(0\M0)\h(p)MPn)}\ 2 . (13) 
Crossing symmetry means 

(0\^ + (0)\h(p),n(p n )) = (n(-p n )\i>+(0)\h(p)) • (14) 

If we make the change p n — ► — p n in the integral (pT3| ) and use (0) we get (remember 
that x — 1/z) 

- z D(z) = q(l/z). (15) 

This relation, which connects the fragmentation function in the physical region < 
z < 1 to the quark distribution in the unphysical region x = 1/z > 1, and vice versa, 
is equivalent to the Drell-Levy-Yan relation Q. 

A further step is needed in order to get the relation (|) . What we have to do is to 
establish a connection between the physical and the unphysical region of the quark 
distribution function (or, equivalently, of the fragmentation function). 

In the physical region < x < 1 the 5-function in (|TTJ) constrains p+ — (1 — x)p + 
to be positive and hence selects positive energy states in the sum over n. Eq. flTT|) 
can be rewritten as 

ff(*) = ^E/ 2 |p0^27T)3 ^ - X P + 'Pn) \(n(Pn)\M0)\h(p)}\ 2 , (16) 
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where p\ = (p^ + Mn) 1 / 2 . The 5-function allows the integration to be simplified giving 

//"CO 
dp n S[p + - xp + -p+] = 2tt / d|p n | \p n \, (17) 
** Pmin 

where 

M 2 (l -x) 2 - M, 

PminW- 2M(1-X) 

We observe now that p m in remains unchanged if we replace x by 2 — x 

Pminix) Pmin(2 x) . (19) 

But Pmin(x) is not the only source of x-dependence in q(x). After exploiting the delta- 
function as in (|17]), the matrix elements appearing in (|TJj) also depend on x. Hence 
Eq. (|16|) is in general non invariant under the substitution x — > 2 — x. However, in the 
large-x limit (which, according to (|18D, is equivalent to the large- \p n \ limit) the matrix 
elements in (|16D tend to become x-independent. The reason is simple. If we describe 
the quarks inside the hadron by Dirac spinors with an upper component u(|p n |) and 
a lower component i>(|p n |), the x-dependence of the matrix elements is contained in 
interference terms of the type u(\p n \)v (\p n \) - see |T^]. Now for large momenta u 



and v must behave as plane waves and their product vanishes when integrated over 
\p n \- Therefore as x gets large the quark distribution tends to become invariant with 
respect to the substitution x — > 2 — x, namely q(x) ~ q(2 — x), or equivalently 

q{l/z)~q{2-l/z), (20) 

in the \axge-z limit. Incidentally, we notice that the same happens in the limit where 
relativistic effects can be neglected. By combining Eq. (|15|) with Eq. (|20"D we finally 
get the relation we have anticipated above 

-D(z) ~ q(2 - l/z). (21) 

Notice that for z > 0.5 this relation connects (approximately) the physical region of 
DIS to the physical region of IA. Eq. (|2l|) is intended to be valid at a fixed and small 
scale /i 2 < 1 GeV 2 . 

We can check the validity of (|20|) by an explicit model calculation. We use a quark- 



diquark model [0 in the framework of the light-cone approach to quark distribution 



functions |16, [17]]. In this model the probability to hit a quark of mass m q and 



transverse momentum k± inside the nucleon, leaving a spectator-diquark with mass 
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mo in the state D, is g_D(x) ~ / d 2 fcj_ \(pjj(x, k 

space wave function of the quark- diquark system with invariant mass M. 



2 , where (Pd{x, k±) is the momentum 

2 _ vr? q +k\ 



+ 



. For the light-cone wave function (Pd{%, k±) we use two different forms: the 



m 2 D +k ± 
1-x 

gaussian type wave function of the Brodsky-Huang-Lepage model [16| and a power- 
law type wave function 



<Pd(x, k±) = A BUL exp(-M 2 /8f3 2 ), 



In Fig. |I] we plot the ratio 



M^k ± ) = A PL (l + M 2 /P 2 ) 
r[z) - 



(22) 



(23) 



(24) 



9(2 - l/z) 

One can notice that r(z) approaches 1 at large z, as we expected, and that the two 
model wavefunctions lead to very similar results. The sharp increase of r(z) as z — > 0.5 
is due to the vanishing of the denominator q(2 — l/z) when its argument tends to zero. 
This is an artifact of the quark-diquark light-cone model which is purely valence-like. 
In more sophisticated models containing a sea of quarks and antiquarks q(x) does 
not vanish as x — > and the increase of r(z) is tamed so that no spurious singularity 
exists z = 0.5. 




Figure 1: The ratio r{z) = q(l/z)/q(2 — l/z) in the light-cone quark model. The solid 
and dashed curves are the results in the light-cone quark model for the Gaussian type 
wavefunction (p2|) and the power-law type wavefunction (^), with m q = 220 MeV, [5 = 450 
MeV, m D = 800 MeV, and a = 3.5. 



Let us come now to the phenomenology of the new relation (21|). If we stick to 
leading order QCD and use Eq. (M), we can translate Eq. fl21|) in terms of structure 
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functions as (remember that at large x or z the evolution is dominated by the quark 
splitting functions) 



F x (z) = 3^(2-1/*); 
zF 2 (z) = -^-L-F 2 (2-l/z) 



(25) 
(26) 



Using standard parametrizations for the DIS structure functions we can predict the 
IA structure functions at large z by means of Eqs. ( p5p6|) . A caveat is in order. Since 
there are only few DIS experimental data for x > 0.7, the quark distributions in this 
region are not very well known. This introduces some uncertainty in our predictions. 

In Fig. we compare the DASP data on z 3 F 2 (z) with the predictions based on the 
new relation ( |26| ) and on the traditional Gribov-Lipatov relation (|5|). For comparison, 
we also show the results for z 3 F 2 (z) based on the approximation (||). For the DIS 
structure functions we used the CTEQ5L parametrization [ I8| . We find that the 
result of the new relation 



is in better agreement with the data at z — > 1. Clearly, 
precision measurements of both F 2 (z) and F 2 (x) at large z and x would allow a more 
conclusive check of (|26|). 
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Figure 2: The structure function z 3 F 2 (z) in e + e~ annihilation. The data are the experi- 
mental results from DASP at Q 2 = 13 GeV 2 |], Rjj. The solid curve is the prediction based 
on ( pT[ ) and (p6|). The dotted curve is the prediction of the traditional Gribov-Lipatov 
relation (^). The dashed curve is the prediction based on (||). For the quark distribution 
functions we used the CTEQ parametrization [18|. 



We already pointed out that Eqs. (^,^) are valid at leading order only. At next- 
to-leading order the evolution of nonsinglet distribution and fragmentation functions 



is different [|P], ^0|. Since the fragmentation functions evolve more rapidly, NLO 
effects lead to a suppression at large z of the LO results for F 2 shown in Fig. 13. Due 
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to the large uncertainty of the present-day data, in this paper we chose for simplicity to 
stick to a leading-order phenomenological treatment. Using NLO splitting functions 
for the fragmentation functions and the new relation ([H]) as the initial condition 
for the evolution one can calculate the NLO corrections to Eqs. (ERES]). 



In conclusion, we presented a new relation between distribution and fragmentation 
functions in their physical regions, which leads to simple testable relations between 
DIS and IA structure functions. A revised form of Gribov-Lipatov relation with a 
color factor and an additional factor of z is also proved to be an approximate relation 
at large z. An immediate application of the new relation connecting q(x) to D(z) 
is in the study of the A polarization near the Z resonance in e + e~ annihilation and 
in polarized lepton DIS scattering. Using the Gribov-Lipatov relation and the QCD 
counting rules for the quark helicity distributions 0], it was found that the data are 
not satisfactorily reproduced Q. We have checked that the situation improves if the 
new relation (plD is used. 
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